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1.  INTRODUCTION 


The  current  stability  prediction  calculations  for  solid  propellant  rocket 
motors  are  based  upon  a  simple  first-order  perturbation  solution  of  the  in- 
viscid  equations  of  motion  for  the  combustion  gas  flow  in  the  chamber.  Thus 
there  is  concern  over  the  applicability  of  the  standard  stability  prediction 
techniques  to  cases  where  high-speed  mean  flows  are  present.  Such  flows  are 
likely  to  occur  in  tactical  rockets  with  large  length- to-diameter  and  low 
port-to-throat  area  ratios.  Nozzleless  motors  represent  an  extreme  example 
of  this  type  of  design;  flow  velocities  may  approach  the  speed  of  sound  at 
the  grain  exit.  High-speed  flow  effects  may  also  be  significant  in  slots 
and  fins  and  around  submerged  nozzles  especially  early  in  the  motor  run  when 
burning  area  is  large  and  flow  channel  area  is  small.  The  purpose  of  this 

t 

study  is  to  assess  the  impact  of  high-speed  mean  flow  on  the  stability  pre¬ 
diction  approach  and  to  generate  correction  techniques  to  be  applied  to  exist¬ 
ing  combustion  instability  models. 

High-speed  mean  flow  affects  the  stability  calculation  in  several  ways. 

The  acoustic  wave  forms  predicted  in  the  linear  model  are  distorted  and  the 
mean  thermodynamic  properties  such  as  the  speed  of  sound  may  vary  significantly 
from  point  to  point  in  the  chamber.  These  variations  may  alter  both  the  pre¬ 
dicted  growth  rates  and  also  the  frequencies  for  the  acoustic  modes.  Correct 
frequency  calculations  are  important  since  they  are  crucial  in  mode  identifi¬ 
cation  in  an  unstable  motor. 

All  acoustic  combustion  instability  models  in  current  use  employ  a  "line¬ 
arization"  of  the  inviscid  equations  of  motion.  This  approach  was  first 
utilized  by  McClure  (Refs.  1  and  2)  and  his  coworkers  and  developed  more  fully 
by  Cullck  (Refs.  3  and  4).  Two  Parameters  representing  respectively  the  ampll- 


1.  Bird,  J.  T.,  McClure,  F.  T. ,  and  Hart,  R.  W.,  "Acoustic  Instability  in  the 
Transverse  Modes  of  Solid  Propellant  Rockets,"  12th  International  Astronautlcal 
Congress,  Academic  Press,  1963. 

2.  Cantrell,  R.  H.,  and  Hart,  R.  W.,  "Interaction  Between  Sound  and  Flow  in 
Acoustic  Cavities:  Mass,  Momentum,  end  Energy  Considerations,"  Journal  of  the 
Acoustical  Society  of  America.  Vol.  36,  No.-  A,  April  196A. 

3.  Culick,  F.  E.  C.,  "Acoustic  Oscillations  In  Solid  Propellant  Rocket 
Chambers,"  Astronautics  Acta.  Vol.  12,  No.  2,  1966. 

A.  Cullck,  F.  E.  C.,  "Interactions  Between  the  Flow  Field,  Combustion,  and 
Wave  Motions  In  Rocket  Motors,"  NWC-TP-53A9,  Naval  Weapons  Center,  China  Lake, 
California,  June  1972. 
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tude  of  the  acoustic  wave  and  the  mean  flow  speed  are  used  In  expanding  the 
pertinent  physical  equations  describing  the  gas  flow.  Combustion  effects  are 
represented  by  admittance  boundary  conditions  at  the  chamber  walls,  nozzle 
entrance  and  head-end  closure.  Only  linear  terms  In  amplitude  and  mean  flow 
Mach  number  are  retained  and  the  result  Is  a  simple  formula  for  the  growth  rate 
which  can  be  evaluated  by  utilizing  linear  acoustic  mode  shapes  in  Integrals 
taken  over  the  surface  and  volume  of  the  chamber  (Ref.  3).  It  Is  clear  that 
this  procedure  is  valid  If  wave  growth  at  the  stability  boundary  only  is  re¬ 
quired,  and  If  an  Incompressible,  uniform  mean  flow  Is  present.  The  assumption 
of  small  wave  amplitude  results  In  what  is  called  an  "tnlttally  valid"  approx¬ 
imation  of  the  growth  process  of  the  waves.  Thus  the  solutions  can  be  consid¬ 
ered  valid  only  at  the  Instant  growth  begins.  It  may  seem  unnecessary  to 
emphasize  this  point  since  prediction  models  indeed  are  based  on  the  stability 
boundary  idea.  However,  an  understanding  is  crucial  If  one  seeks  to  Improve 
the  linear  model  by  evaluating  higher-order  terms  In  the  asymptotic  expansions. 
Specifically,  the  linear  model  assumes  that  the  small  parameters  are  related  by 
the  limit  process 
1  Im  (e/M.  )  ■  0 
e,  Mfa-  o 

where  c  is  the  wave  amplitude  (proportional  to  the  Mach  number  of  the  acoustic 
velocity)  and  Mfa  Is  the  mean  flow  Mach  number  at  the  burning  surface.  It  Is 
readily  demonstrated  that  if  one  desires  to  assess  the  effects  of  mean  flow 
compressibility  and  spatial  variation  of  the  mean  flow  field  thermodynamic 
properties  then  higher-order  terms  (at  least  to  0(Mb  ))  must  be  retained  despite 
the  fact  that  Mb  at  the  surface  Is  very  small.*  In  order  that  the  expansion 
process  does  not  lead  to  spurious  results,  it  Is  then  necessary  to  modify  the 
limit  process  given  In  Equation  1  to  the  following: 

1 im  (c/Mfa2)  -  0  {2) 


That  Is,  c  ■  o(Mb  )  In  the  notation  of  perturbation  theory.**  In  view  of 
*  Mean  flow  Mach  number  at  the  burning  surface  Is  typically  In  the  range 
0.001  <  M.  <  0.01. 

D 

**  If  this  restriction  is  not  made,  then  it  Is  necessary  to  evaluate  terms  of 
2 

order  c  H.  to  produce  a  valid  asymptotic  expansion. 

D 


the  small  magnitude  of  M^,  this  represents  a  severe  restriction  on  the  ampli¬ 
tude  of  the  wave.  The  difficulty  appears  when  one  attempts  to  correlate 
experimental  results  by  use  of  this  theory.  By  definition,  the  data  repre- 

2 

sent  finite  waves,  and  even  the  smallest  dtscernable  wave  motions  imply  o>Mb 
which  violates  the  fundamental  assumption.  Limitations  of  the  sort  described 
can  only  be  overcome  by  either  reverting  to  a  completely  numerical  solution 
of  the  problem  or  by  use  of  a  more  appropriate  perturbation  scheme.  In  regard 
to  the  latter  approach,  It  appears  that  a  modification  of  the  two-variable  method 
or  method  of  multiple  scales  will  eventually  allow  generation  of  analytical 
solutions  which  are  valid  for  a  less  limited  range  of  values  of  wave  amplitude 
(Refs.  5  and  6).  Such  expansions  can  yield  uniformly  valid  approximations 
which  would  be  of  great  utility  In  understanding  the  self-1 Imltlng  character¬ 
istics  or  "limiting  amplitude"  behavior  of  solid  rocket  instability.  Such  an 
approach  was  foregone  in  the  present  study,  since  emphasis  was  to  be  on 
devising  computer  codes  which  are  now  in  widespread  use.  Thus  the  expansion 
process  implied  by  Equation  2  forms  the  basis  for  the  analysis  given  here. 

All  dynamic  and  thermodynamic  parameters  of  the  mean  and  oscillating  flow 
fields  are  expanded  In  both  c  and  K^.  Figure  1  will  aid  the  reader  In  Inter¬ 
preting  the  degree  to  which  we  are  attempting  to  Improve  the  linear  theory. 

Retention  of  terms  of  order  (e,  M.)  yields  the  linear  stability  model.  In 

®  2 

what  follows,  we  extend  the  model  to  Include  effects  of  0(c,  Mb  )  In  order  to 
assess  the  influence  of  nonuniform,  compressible  mean  flow  and  mean  chamber 
thermodynamics.  As  in  the  "classical"  model,  combustion  effects  and  the 
Influence  of  the  nozzle  are  represented  by  appropriate  admittance  boundary 
conditions.  Also,  viscosity  is  assumed  to  be  negligible  except  Insofar  as  it 
affects  the  wave  motions  near  the  solid  boundaries  of  the  chamber.  These 
influences  can  be  accounted  for  by  corrections  to  the  surface  admittance 
(Ref.  7).  However,  vortlclty  generated  at  the  boundary  may  be  convected  through¬ 
out  the  flow  field  and  subsequently  affect  both  wave  growth  and  oscillation 
frequency  (Refs.  8  and  9).  Since  these  effects  may  be  especially  Important 
in  the  presence  of  high  mean  flow  Mach  numbers,  their  influence  Is  retained  in 

5.  Cole,  J.  D.,  and  Kevorkian,  J.,  "Uniformly  Valid  Asymptotic  Approximations 
for  Certain  Nonlinear  Differential  Equations."  Nonlinear  Differential  Equations 
and  Nonlinear  Mechanics,  Academic  Press,  1983. 

T.  Nay f eh.  A.  H.,  Perturbation  Methods,  Wiley,  1973. 

7.  Flandro,  G.  A.,  "Solid  Propellant  Admittance  Corrections,"  Journal  of  Sound 
and  Vibration,  Vol.  38*  Ho.  3,  197A. 


the  formulation.  The  analysis  Is  carried  out  for  the  full  three-dimensional 
geometry  with  arbitrary  shape.  The  results  are  applied  to  longitudinal  waves 
In  a  high  L/D  cylindrical  grain  to  demonstrate  their  potential  relevance  in 
the  stability  prediction  problems.  It  Is  shown  that  significant  changes  in 
mode  frequency  and  growth  rates  occur;  the  corrections  grow  quadratlcal ly 
with  increases  in  the  port  length-to-diameter  ratio  L/D. 


8.  Culick,  F.  E.  C.,  "Rotational  Axlsymmetrlc  Mean  Flow  and  Damping  of 
Acoustic  Waves  In  a  Solid  Propellant  Rocket,"  AIAA  Journal ,  Vol .  k.  Mo.  8, 
August  1966. 

9.  Flandro,  G.  A.,  "Rotating  Flows  In  Acoustically  Unstable  Rocket  Motors," 
PhD  Thesis,  California  Institute  of  Technology,  1967. 
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2. 


ANALYSIS 


In  what  follows,  the  equations  of  motion,  boundary  conditions,  and 
assumptions  which  form  the  basis  for  the  analysis  are  first  reviewed.  The 
mathematical  strategy  is  then  formulated  and  applied  to  the  problem.  Results 
are  given  in  completely  general  form.  The  reader  more  interested  in  the 
applications  than  the  details  of  the  calculation  might  skip  all  of  this 
section  except  the  review  of  the  assumptions. 


2.1.  Assumptions 

A  rigid  internally  burning  propellant  grain  is  assumed.  Only  regions 
of  the  flow  which  are  subsonic  are  considered.  Influence  of  the  sonic  nozzle 
throat  are  represented  by  appropriate  admittance  functions.  Although  an 
inviscid  fluid  is  assumed  (no  viscous  force  terms  are  retained  in  the  momen¬ 
tum  equations)  the  flow  is  allowed  to  be  rotational;  thus  a  realistic  mean 
flow  field  is  accommodated  (Ref.  8)  and  vorticlty  generation  and  trans¬ 
port  is  represented.  The  combustion  process  is  assumed  to  be  concentrated 
at  the  bounding  surface  of  the  chamber  and  its  sensitivity  to  acoustic  fluc¬ 
tuations  is  represented  by  an  admittance  boundary  condition  as  in  the  linear 
stability  analysis.  The  effect  of  regression  of  the  burning  surface  as 
propellant  is  consumed  is  neglected  since  the  associated  characteristic  time 
is  long  compared  to  the  period  typical  of  the  gas  oscillations  and  amplifi¬ 
cation  time  for  the  waves.  The  gas  is  assumed  thermodynamically  perfect. 

2.2.  Equations  of  Motion 

.  The  motion  of  a  compressible  inviscid  gas  is  governed  by 


Iy  +  V  •  (pu )  =  o 


(3) 


*  Some  effects  of  rigid  body  motion  of  the  combustion  chamber  on  structure 
and  growth  of  acoustic  waves  are  assessed  in  Reference  9- 


8 


D-  VP  -  o  (4) 

p  Dt  +~ 


P 


21 

ot 


DP 

Dt 


o 


(5) 


P  -  pT 


(6) 


The  equations  are  written  in  terms  of  the  dimensionless  variables: 


P  -  P'/Po 

p  =  p'/po 


u  =  u'/a 
—  —  o 

T  -  T'/To 
r  -  £'/R 
t  =  (ao/R)f 

where  R  is  a  characteristic  chamber  radius  and  aQ  is  the  stagnation  speed  of 
sound  in  the  chamber.  Dimensional  quantities  are  denoted  by  primes  and  sub¬ 
script  o  indicates  the  stagnation  values  of  the  principal  thermodynamic  vari¬ 
ables  (dimensional)  in  the  absence  of  wave  motion. 

2.3.  Perturbation  Expansions 

In  order  to  extend  the  linear  theory  it  is  necessary  to  utilize  double 
perturbation  series  of  the  principal  variables  in  terms  of  the  two  fundamental 
scaling  parameters  M^  and  e.  M^  is  the  mean  flow  Mach  number  at  the  burning 
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**••■»*  ixr.aii’* 


\ 

l  t 


p^o)  =  !  +  Mfa2  P  +  0  (M^) 


p(l)  =  p(IO)  +  M  pO  1)  +  0  (M  2) 

D  D 


-u<10>  ♦  .0  («h2) 
—  —  D —  D 


P  reflects  the  effects  of  compressibility  on  the  steady  pressure  distribution 
in  the  cavity.  P  and  U  are  usually  determined  as  functions  of  position  in 
the  chamber  as  part  of  the  motor  internal  ballistics  design  process.  For 
complicated  geometries  they  cannot  be  described  in  terms  of  simple  mathemati¬ 
cal  functions,  but  can  be  handled  numerically  in  the  calculations. 

2.  A.  Boundary  Conditions 

Since  the  combustion  process  takes  place  in  a  thin  lamina  of  fluid  at  the 
wall  of  the  chamber,  it  is  possible  to  represent  the  effects  of  combustion  as 
a  boundary  condition  on  the  flow  at  the  interface.  This  is  accomplished  by 
introducing  the  admittance  function  A  such  that 


n  •  u  =  -M.  (— )A 

-  by 


at  the  boundary,  ui  and  P  are  the  (complex)  velocity  and  pressure  fluctuations 
at  the  wall,  n  is  an  outward  pointing  unit  vector  normal  to  the  chamber  wall, 
and  y  the  ratio  of  specific  heats  is  inserted  for  convenience.  This  notation 
for  A  is  the  same  as  originally  employed  by  Culick  (Ref.  3).  In  later  work 
he  altered  this  definition  somewhat  by  defining  n  to  be  an  i nwa rd  pointing 
normal  unit  vector  and  by  omitting  the  scaling  factor  from  the  definition. 
To  distinguish  the  two  definitions,  we  put 


"b 


1 


-u'-  ^  ■ 


where  is  the  admittance  function  normally  used  in  the  literature  (Refs.  4 
and  10).  The  use  of  the  older  notation  is  adhered  to  in  this  analysis  since 
higher-order  expansions  in  are  evaluated--carrying  superfluous  M^'s  in 
the  already  sufficiently  complicated  algebraic  expressions  is  clearly  not  a 
sensible  procedure.  It  is  well-known  that  the  algebra  involved  in  perturba¬ 
tion  expansions  increases  enormously  as  higher  order  terms  are  retained;  it 
is  thus  essential  to  use  efficient  notation.  An  additional  advantage  of 
the  present  notation  is  that  numerical  values  of  | A |  are  in  the  range 
- 1 0<  j  A  j  <10  such  that  very  tiny  numbers  and  resulting  additional  possibility 
of  error  does  not  affect  the  calculations.  Note  that  A  is  a  complex  number: 

A  =  A^  +  l  A^  (13) 


It  will  be  shown  that  knowledge  of  both  real  and  imaginary  parts  is  necessary 
in  the  proper  evaluation  of  high-speed  mean  flow  effects. 

2.5.  Derivation  of  Wave  Equation 

Equations  (3)  through  (6)  may  be  recombined  to  yield  the  set 


DP 


3t 


+  yPV-U  =  U-VP 


(14) 


VP 


p  df  +  r =  0 


(15) 


Utilizing  expansions  (9)  and  (10)  one  finds  to  0(e)  (retaining  terms  to  0  (M,Z)) 
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|£  +YV.u(1)  =  -Mb  [U-Vp(,)]  -  Mb2  [u(,)-Vp+  YpVu{,)]  (16) 


10.  Anonymous,  T- BURNER  MANUAL.  CPIA  Publication  No.  191,  November  1969- 
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MK 
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Y  Y 


+  P(,)v»p] 
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where  ^  =  7  X  is  the  mean  flow  vorticity. 

The  wave  equation  governing  the  acoustic  pressure  fluctuations  is  derived 
oy  subtracting  the  divergence  of  the  momentum  equation  (17)  from  the  time 
derivative  of  the  continuity  equation  (16): 


>2„(1) 


-V2P(,)  =  M  t-|r  (U*Vp(,))  +  YV*V(U*u(,))  -  Y7*(u(,)  X  C)J 
at  d  at  --  ~ 


(18) 


‘  Mb2  llr  (u(1)^p+rpv-u(1))  +  ^•(EJ&(1)  +  P(1)vp)l 


b  lSt 


,  2 

Again,  terms  of  order  are  retained  to  correctly  represent  the  influence  of 
the  high-speed  compressible  mean  flow.  It  is  important  to  note  that  this 
implies  that  all  0  (M^2)  corrections  (e.g.  P^^,  ^ ,  etc.)  must  eventually 

be  calculated  to  properly  evaluate  the  wave  equation.  In  anticipation  of  oscil¬ 
lations,  it  is  appropriate  to  assume  that  all  dependent  variables  exhibit  expo¬ 
nential  time  dependence: 
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Y  Is  Inserted  In  equation  (19)  for  convenience;  <j^  is  spatial  distribution 
of  acoustic  velocity,  and  K  is  the  complex  frequency 

K  =  Si  +  -tA  (21 ) 


13 


14 


boundary  condition.  It  is  useful  at  this  point  to  further  expand  this  set  in 


terms  of  the  mean  flow  Mach  number  M^.  Note  that  the  eigenvalue  K  must  also 

2 

be  expanded;  the  goal  of  the  calculation  is  to  find  K  correct  to  0(M.  )  to 

D 

represent  growth  rate  and  frequency  shifts  due  to  the  high-speed  mean  flow. 

Put 


p('>.p<,‘».M(,p<">.H2p"2>  .... 
D  D 


(26) 


.„„(„<">  .a(">)  .K,2  to"2)  .*a<I2>)  .... 


9(,)  -,<><»  .. 


h<"  .(,<»»  .Mhh<">  .... 
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p^10^  is  the  unperturbed  mode  shape  of  the  acoustic  wave;  p^1^  etc.  represent 
distortions  of  the  waveform  caused  by  the  mean  flow  perturbations.  Inserting 
expansions  (26)  into  the  wave  equation  (22)  and  boundary  condition  (24) ,  one 
finds  a  series  of  boundary  value  problems  in  increasing  orders  of  M^;  the  0  (1) 
problem  obviously  corresponds  to  unperturbed  acoustic  wave  motion  in  the 
chamber.  Thus 


'Vp(l0>  .„(">»  p<"» 


(27) 


n*Vp^®^  *  0  on  bounding  surfaces 
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Previous  comments  on  increasing  algebraic  complexity  with  order  of  expansion 
are  clearly  appropriate.  Expressions  for  the  velocity  amplitude  (and  its 

expansion  in  Mach  number  M^)  are  found  in  terms  of  pressure  fluctuation  p^  by 
use  of  the  momentum  equation.  Thus 


(10)  _  -tVp^10) 


and  3(n)  -rnoT 


E(">  -  (°("(io) 'a(">)  vp(,0) 


+  7noT  <vp<'0)  Xi) 


2.6.  Zeroth-Order  Calculations 

Equations  (27)  and  (28)  represent  the  classical  acoustics  of  the  combustion 
chamber.  E/act  solutions  can  be  found  for  simple  chamber  geometry;  numerical 
techniques  provide  information  on  mode  shape  p^^  and  frequency  for  the 

sometimes  complicated  configurations  used  in  solid  propellant  rockets.  Since 
this  phase  of  the  problem  is  thoroughly  documented  in  the  literature,  no  more 
need  be  said  at  this  point  except  that  the  validity  of  the  higher-order  growth 
rate  and  frequency  calculation  obviously  depends  on  having  correct  estimates  for 
p^10^  and  for  each  chamber  mode.  The  unperturbed  frequency  is  found  from 

the  eigenvalue  in  dimensional  form  (cycles/unit  time)  to  be 


(38) 
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2.7.  First-Order  Calculations 

Equations  (29)  and  (30)  govern  what  is  commonly  referred  to  as  the  "linear 
stability"  model  of  combustion  instability.  In  this  calculation,  one  is  seldom 
interested  in  the  0 (M,  )  waveform  distortion  p^'^;  what  is  required  is  an  esti- 
mate  for  the  growth  rate  A'  1  and  frequency  correction  .  Culick  (Ref.  3) 

derived  expressions  for  these  quantities  by  the  Green's  function  method.  A 
more  direct  method  is  to  multiply  (27)  by  P^''^  and  subtract  it  from  (29) 
multiplied  by  the  unperturbed  mode  shape  p^*^.  Integrating  the  result  over 
the  chamber  volume,  employing  the  divergence  theorem  and  the  boundary  con¬ 
ditions  (28)  and  (30)  yields  after  a  few  lines  of  algebra: 


(fl(,)  +  cA(,)) 


2fi^°^E2 


g(,CV,0)dV  + 


h<10>P<,0>dS 


(39) 


where  g^^  and  h^^  are  given  by  equations  (33)  and  (3k),  and  the  normaliza¬ 
tion  constant  E  is  found  from 


E2  = 


p(,0>  dV 


(*•0) 


II 


Equation  (39)  can  be  further  simplified  (see  the  papers  by  Culick  for  details) 
by  application  of  vector  integral  theorems  with  the  results: 


Growth  rate  (dimensional)  a  is  found  from  (42)  to  be 


a 


(1)  .  . 
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(43) 


where  R  is  the  scaling  length  for  the  chamber;  aQ  is  the  speed  of  sound. 

Equation  (43)  forms  the  basis  for  all  stability  prediction  models.  The  results 
are  discussed  at  length  in  the  literature  (c.f.  Refs.  3,  4,  8,  10,  etc.)  and 
no  more  need  be  said  except  that,  again,  it  is  obvious  that  higher-order  ap¬ 
proximations  depend  on  correct  estimates  for  A ^  and  0^.  The  0  (M.  )  fre- 

(1 )  D  . 

quency  correction  Q  is  usually  ignored  in  stability  calculations.  This  is 

probably  due  more  to  lack  of  information  regarding  the  imaginary  part  of  the 

admittance  than  to  it  being  a  negligible  contribution.  The  frequency  correct ioi 

in  cycles/unit  time  is  found  from  to  be 
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■<Kii>  Mb  Q 
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(44) 


*  The  superscript  is  used  to  emphasize  the  order  of  the  estimate;  higher-order 
approximations  are  the  goal  of  this  study. 


This  expression  is  evaluated  in  the  results  section  of  the  report  for  typical 
sets  of  parameters  to  establish  its  potential  importance. 


Another  part  of  the  first-order  formulation  is  required  before  the  0  (M,z) 

2  b 

corrections  can  be  determined.  Examination  of  the  M.  equation  shows  that  one 

b 

must  know  the  mode  shape  correct  to  0  (M,  )  before  the  higher-order  equations 

“  (|i) 

can  be  evaluated.  That  is  the  waveform  correction  p  must  first  be  deter¬ 
mined.  There  are  two  approaches  to  this  problem.  One  is  to  seek  the  eigen¬ 
functions  for  equation  (29)  directly.  This  can  be  done  analytically  only  for 
simple  geometries.  This  method  is  demonstrated  in  a  later  section  of  the 
report.  A  more  general  method  is  to  use  a  Green's  function  solution  of 
equation  (29)  using  expansions  in  terms  of  the  unperturbed  eigenfunction  (Refs. 


3  and  4).  This  yields  an  infinite  series  representation  for  p 
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g(,0)p(10)dV  + 


h(,0V’0)dS 


(45) 


where  N  refers  to  the  particular  mode  of  oscillation  in  question;  6  identifies 
each  of  the  remaining  infinity  of  modal  eigenfunctions.*  It  is  the  experience 
of  the  writer  that  this  is  usually  a  very  slowly  convergent  series.  It  is 
possible  that  direct  numerical  solutions  of  equation  (29)  will  be  a  more  ef¬ 
ficient  technique  for  determining  the  mode  distortion  effects. 

2.8.  Second-Order  Corrections 

We  finally  arrive  at  the  target  of  the  present  study--to  determine  the 

2 

effects  of  mean  flow  to  0  (Mfa  )  on  stability  and  frequency.  It  is  important 
to  note  that  no  influence  of  mean  flow  compressibility  or  variation  of  mean 

*  For  three  dimensional  chambers  N  (and  3)  are  three- i nteger  sets  which  identify 
the  mode. 
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chamber  thermodynamic  parameters  ts  reflected  In  the  results  of  the  first 
order  calculation;  no  questions  regarding  high-speed  mean  flow  effects  can 
be  answered  by  use  of  Equations  (4l)  and  (42).  We  now  calculate  ft^  and 
which  reflect  the  effects  of  Interest.  Multiplying  Equation  (27)  by  p^^ 
and  subtracting  (31)  multiplied  by  p^®^  yields 
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after  Integration  over  the  chamber  volume  and  application  of  the  boundary  con¬ 
ditions.  Evaluation  of  these  expressions  gives  the  required  Information.  It 

(12) 

Is  not  necessary  to  determine  p  explicitly  unless  even  higher  order  cor¬ 
rections  are  required. 

The  application  of  these  results  is  demonstrated  in  the  following  section 
for  a  typical  motor  configuration.  It  Is  Important  to  notice  that  the  form  of 
the  frequency  and  growth  rate  corrections  Is  identical  to  the  first  order  ex¬ 
pression.  The  significance  of  this  observation  is  that  the  very  same  type  of 
calculation  as  used  in  the  linear  stability  model  is  required  to  determine  the 
corrections.  Once  the  waveform  distortion  Is  known,  the  rest  of  the  calcula¬ 
tion  Is  handled  by  simply  adding  correction  terms  to  the  Integrals  already 
evaluated  in  the  standard  codes.  The  expressions  are  somewhat  more  complicated 
algebraically  than  their  first-order  counterparts,  but  they  are  easily  Incor¬ 
porated  by  simply  adding  them  to  existing  codes.  In  addition  to  the  unperturbed 
mode  shape  data  utilized  In  the  standard  codes,  one  must  provide  the  following 
additional  input: 

1.  Waveform  distortion  (p^^) 

2 

2.  0  (M^  )  pressure  distribution  (?) 

3.  Reasonable  estimate  of  mean  flow  velocity  distribution  (U) 

2 

The  appendix  Includes  a  discussion  of  methods  for  estimating  the  OfM^  ) 
compressible  pressure  distribution  and  mean  flow  velocity  pattern  In  a  tubular 
grain. 
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3.  APPLICATIONS 


The  purpose  of  this  section  is  to  demonstrate  in  detail  the  application 
of  the  theory  set  forth  in  the  previous  parts  of  the  report.  A  simple  motor 
geometry  is  analyzed  to  enable  analytic  determination  of  all  quantities  of 
interest. 

The  simplest  configuration  which  exhibits  the  features  we  desire  to 
emphasize  is  a  tubular  grain  with  high  L/D  (length-to-diameter)  ratio.  In 
what  follows,  the  effects  of  high-speed,  compressible  mean  flow  on  the  fre¬ 
quency  and  growth  rate  of  longitudinal  waves  in  such  a  motor  are  determined; 
the  results  are  directly  applicable  to  several  tactical  rocket  designs  cur¬ 
rently  under  development.  Figure  2  shows  the  geometry  assumed  and  the  co¬ 
ordinate  system.  The  Appendix  contains  a  detailed  analysis  of  the  mean  flow 
pattern  and  thermodynamics  appropriate  to  this  configuration. 

Assuming  longitudinal  modes  of  oscillation  (3  =  l  =  I,  2,  3-..),  solution 
of  the  zeroth-order  equations  yields  the  familiar  results: 
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where  l  is  the  mode  integer  (l  =  1  is  the  first  longitudinal  mode,  etc.).  The 
unperturbed  frequency  (from  eqn.  38)  is 
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Converting  to  the  more  familiar  admittance  notation  we  find: 
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(^)(Ab(r)  -  Hb)  +0  (Hb2) 


(57) 


or 


a 


<*)<V2Mb)  +  0  <Mb2> 


(58) 


where  Rb  =  A^1^  +  Mb  is  the  response  function  (Ref.  10).  In  an  actual  motor 
calculation  these  results  would  be  supplemented  by  nozzle  losses,  particulate 
losses,  etc. 
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To  proceed  further,  it  is  necessary  to  determine  the  waveform  distortion 
by  solving  equation  (29)  for  p^'^.  For  the  cylindrical  geometry,  the  first- 
order  velocity  distribution  is 


-  i.  s  i  n 
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where  e  is  the  axial  unit  vector  and  function  g 
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Thus  the  differential  equation  for  p 


d!_P.(,,)  +  k  2  (11) 

2  +  k(i  P 

dz  4 


=  -^-i-k,, 


(11) 


reduces  to: 
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and  the  boundary  conditions  (from  equation  (30))  reduce  to: 
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(Nozzle  entrance) 


Solutions  are  found  easily  by  conventional  methods  ana  trie  result 
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figure  3  shows  plots  of  the  waveform  p  =  p  '  -  1,  p  '  '  for  t  ie  'irst  wo 

(-o'1  b 

long'  tuc!  ;na’  modes.  A  ,  the  imaginary  part  o'  the  surface  aoml  rvner,  ins 

only  a  small  effect  on  the  real  part  of  the  mode  ’.nop."  the  p  vs  w-  non  a 

w  tr  A'  =  1,  a  typical  value,  and  severa  vi'ie-.  n;  *ho  r>  a  «*  trv  nd 

(r' 

m.  ttance,  A  ‘ .  -rhe  most  :mportant  feature  of  the  •'esu’ts  is  toe  generrt  on 
of  the  strong  0  (M^)  maginary  part  of  the  moce  shape.  As  w.’’  be  shown  this 
exeats  n  strong  in'’uence  on  the  chamber  f  recurrr.  No*o  thr.  there  i  <  «- 

•  -«!i  i '  leant  sens'  t'v:ty  o'  the  imaginary  mod  d ' «  o:  or  *.<>  thr  1  pr  i  i  o' 
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Function  g  for  the  simple  geometry  assumed  is 
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*ii  algebra.  in  tern’s  of  the  same  notation,  function  n*  i 


i  <1  -  *«">)  r  ,*<">  . 


It '  '  1 ;  =”tC  (A vr  ^  r  cos(k.z)  -  < .  ( A  ^ 1  +  tA^F)B  k.z  sin  (k  z) 


(r)  ,(<),,  2  2  „  , 
t  t  t.t\  :  k  Z  COStk,/) 


j  2  s  i  fi  f  < ^  z )  +  ii?4")  +  k  '  3;  KjZ  cos(kv/) 


+  [6  -  ■/!  ■  2 <  B j  k„2  z‘  sin(k.z) 


+  2  k  5  z-5  cos  (k.z) 


•n  boundary  surfaces. 


The  integrals  can  now  be  performed  and  (after  some  tedious  algebra!)  one 
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where  the  previously  calculated  values  of  ^ ^  and  '  have  been 
I  i  zed .  Presence  of  dampinq  would  thus  alter  these  results,  but  they  ere  r 
sentative  of  the  influence  of  the  high-speed  mean  flow.  Converting  to  dim. 
sional  form,  the  growth  rate  and  frequency  corrections  are 


t  i  nos 
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Growth  rate  and  frequency  are  corrected  by  adding  (70)  and  (71)  to  (55)  and 
(56)  respectively. 

Plots  of  these  expressions  are  displayed  in  Figures  **  through  8  for  first, 
second  and  third  longitudinal  modes.  Note  that  the  growth  rate  correction  is 
linearly  dependent  on  the  imaginary  part  of  the  wall  admittance  function  A ^ . 
Since  A^  is  positive  for  lower  frequencies  typical  of  longitudinal  oscilla¬ 
tions,  this  is  an  additional  source  of  acoustic  gain;  for  higher  frequencies, 
high-speed  flow  contributes  additional  damping.  The  magnitude  of  the  growth 
rate  increment  is  also  dependent  on  the  real  part  of  the  admittance  as  illus¬ 
trated  and  is  proportional  to  the  square  of  the  port  L/D  ratio.  It  is  clear 
that  this  could  be  an  important  growth  rate  element  in  motors  with  high  L/D. 

For  instance,  in  a  tactical  motor  with  L/D  =  20,  M.  =  0.007,  Y  ■  1.2,  A^  =  2, 
,  • ,  b 

A  c  =  2  and  f  *  300  Hz,  the  growth  rate  estimate  would  be  increased  by  about 
^5  sec  '.  As  shown  in  Figures  6  through  8,  the  frequency  shift  depends  on 
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and  the  square  of  the  motor  L/D.  A  tactical  motor  with  the  parameters 
just  mentioned  would  experience  a  frequency  shift  of  about  -90  Hz  (30  percent 
change)  which  is  clearly  significant.  Since  effects  of  damping  phenomena  have 
not  been  included  in  these  numerical  calculations,  the  growth  rate  modification 
and  frequency  shifts  are  most  likely  not  this  large  in  the  actual  motor. 


Figure  4.  Growth  Rate  Increment-First  Longitudinal  Mode 


4. 


CONCLUSIONS 


The  analysis  presented  in  this  report  demonstrates  the  importance  of 
several  factors  which  are  not  included  in  standard  solid  propellant  rocket 
combustion  stability  predictions.  The  formal  linearized  model  of  acoustic 
instability  is  based  op  he  i«e,ur.:>t  ion  thu  the  moan  p'essure  f  i  e  >  ri  is  uniforr: 
and  that  the  I  low  is  incompressible.  [he  present  analysis  shows  that  in  high 
L/D  rockets  or  in  any  motor  configuration  which  exhibit  high  speed  mean  flow, 
tfv  compressibility  and  mode  d  i  >r  t  i  on  of  t  s  rt’osi  bo  accounted  ''or  .  Fre¬ 
quency  shifts  of  as  much  as  3U  percent  of  the  linear  acoustic  frequency  p re- 
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APPENDIX 


COMPRESSIBLE,  ROTATIONAL  AX  I  SYMMETRIC 
FLOW  IN  A  SOLID  PROPELLANT  ROCKET 


1.  INTRODUCTION 


Required  for  the  assessment  of  high-speed  flow  effects  in  acoustic  com¬ 
bustion  instability  is  a  practical  and  realistic  model  of  the  compressible  mean 
flow  in  the  combustion  chamber.  Such  a  model  would  also  be  of  obvious  utility 
in  studies  of  erosive  burning,  growth  of  metallic  oxide  particles  and  velocity 
coupling  effects.  Results  from  detailed  computer  models  of  the  chamber  flow 
are  usually  not  in  a  form  useful  in  theoretical  investigations;  one-dimensional 
models  do  not  represent  the  geometrical  flow  features  which  are  crucial  in  the 
problems  mentioned. 

The  Incompressible  axisymmetrlc  solution  for  steady  rotational  flow  in  a 
cylindrical  rocket  grain  devised  by  Culick  (Ref.  11)  is  the  starting  point 
for  the  present  investigation.  It  was  shown  by  Dunlap  et.  al.  (Ref.  12)  that 
the  Culick  model  closely  satisfies  the  viscous  equations  of  motion,  and  the 
results  were  experimentally  verified  using  a  cold  flow  apparatus.  In  what 
follows,  Culick's  solution  is  generalized  to  Include  the  effects  of  compressi¬ 
bility  In  the^mean  flow  field.  The  results  of  the  extended  analysis  are  veri¬ 
fied  by  comparison  to  a  one-dimensional  numerical  solution.  An  Important  feature 
of  the  model  is  its  simplicity;  this  greatly  enhances  its  value  in  theoretical 
studies  of  the  type  described  above. 

£riLo&Ui,«G  £>.uK  NOT  FIiMlJO 

2.  ANALYSIS 

The  steady,  compressible,  rotational  axially  symmetric  flow  of  a  perfect 
gas  is  governed  by 
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Propellant  Rocket  Motor,"  AIAA  Journal.  Vol.  12,  No.  10,  Oct.  197*».  pp.  1AA0-1AA2. 
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where  the  continuity  equation  has  been  exactly  satisfied  by  specification  of 
the  stream  function  p  such  that 
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u  and  w  are  the  radial  and  axial  velocity  components  respectively.  Equation 
A1  is  the  vorticity  equation;  c  is  the  azimuthal  component  of  vorticity.  Radial 
and  axial  vorticity  components  are  zero  in  axi symmetric  flow.  Equation  A2  is 
the  momentum  equation  and  Equation  A3  is  the  energy/state  relationship.  Viscous 
force  terms  do  not  appear  explicitly  in  the  momentum  equation  because  the  ap¬ 
propriate  length  scale  is  the  chamber  radius.  This  makes  the  viscous  force 
much  smaller  than  momentum  and  pressure  effects;  no  boundary  layer  of  the 
classical  type  is  formed  and  vorticity  fills  the  entire  flow  field.  All  flow 
properties  are  made  dimensionless  in  terms  of  the  stagnation  thermodynamic 
properties  a^,  Po  and  po  at  the  head-end  of  the  chamber;  chamber  radius  R  is 
the  characteristic  length.  Boundary  conditions  are 
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The  latter  constraint  requires  the  flow  to  satisfy  the  no  slip  condition  at 

the  burning  surface;  this  introduces  vorticity  into  the  flow  and  also  forces 

2 

it  to  behave  as  a  viscous  fluid  . 

The  problem  is  readily  solved  by  standard  perturbation  techniques;  the 

I 

mean  flow  Hach  number  at  the  burning  surface  is  the  appropriate  small 
parameter.  The  expansions  are 
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The  first  order  solution  is  governed  by 
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where  to  satisfy  the  boundary  conditions  specified.  This  differs 

from  the  vorticity  used  by  Culick  (Ref.  11)  by  an  additional  factor  of  ir.  The 
solution  to  Equation  A8  is 
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and  the  corresponding  velocity  components  are 
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Note  that  the  next  terms  in  the  expansion  for  are  0  (M  so  that  Equation 
A10  represents  the  velocity  field  to  a  high  degree  of  precision. 

Equation  A2  is  now  expanded  to  yield  information  regarding  the  thermody¬ 
namic  properties  of  the  flow  field. 

2 

Retaining  terms  of  0  (M^  ),  one  finds 
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Noting  that  ip^°^ 
wi th  the  resul t 


/equation  All  can  be  integrated  directly 
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The  constant  of  integration  is  zero  to  satisfy  stagnation  requirements  at  the 
head-end.  Inserting  Equation  A9  yields 


and  using  Equation  A13  in  an  expansion  of  Equation  A3  gives 


These  functions  represent  models  of  the  compressible  thermodynamic  prop- 

1, 

erties  correct  to  0  (M^  ). 

3.  RESULTS 

The  analysis  yields  a  simple  and  practical  representation  for  the 
flow  field  in  a  tubular  rocket  grain. 

Rewriting  the  equations  in  terms  of  physical  variables, 
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where  v.  *  p  r/p  -a  M,  is  the  velocity  at  the  chamber  boundary  due 

D  p  O  O  D 

to  efflux  of  combustion  gases,  r  and  z  are  the  radial  and  axial  position 
c  coordinates  and  R  Is  the  chamber  radius. 


Figure  A1  shows  the  streamline  pattern  corresponding  to  the  above  so¬ 
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Normalized  axial  position 


Figure  A2  illustrates  the  velocity  profiles.  The  upper  plot  shows  the 
axial  velocity  distribution;  the  lower  shows  the  radial  velocity  profile.  The 
axial  velocity  is  conveniently  referenced  to  the  centerline  speed 
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while  the  appropriate  reference  speed  for  the  radial  component  is  the  efflux 
veloci ty  v^. 

Figure  A3  shows  the  pressure  variation  along  the  motor  axis.  It  is 
apparent  from  Equations  A15  that  the  pressure  varies  quadratical ly  with 
axial  position.  At  the  centerline,  the  variation  is  given  by 


Pc-Po 
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where  the  pressure  scale  is  the  reference  dynamic  pressure  based  on  chamber 

pressure  at  the  head-end  and  the  burning  surface  Mach  number  M.  .  A  typical 

2  2  ° 
value  for  yPoMb  /2  is  0.02  (corresponding  to  Po  »  600  lb/in  and  M^  ■  0.007). 

At  the  burning  surface,  the  pressure  variation  is 


(A  18) 


Thus  the  radial  pressure  gradient  is  much  smaller  than  the  axial  component. 


Since  the  flow  field  predictions  based  on  the  present  analysis  are  not 
valid  when  sonic  speed  is  reached.  It  is  necessary  to  calculate  the  limiting 
axial  position.  It  is  readily  shown  that  a  Mach  number  of  unity  is  exceeded 
if  z/R  exceeds  the  critical  value 
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Chamber  axis 


where  M^  is  the  burning  surface  Mach  number  and  y  is  the  ratio  of  specific  heats. 
For  example,  if  y  =  1.2,  M^  =  0.007,  (z/lt)  *  43.4. 

Thus  for  these  typical  conditions,  sonic  flow  is  reached  if  the  motor  L/D  ex¬ 
ceeds  about  22. 

4.  HIGHER-ORDER  CORRECTIONS 


The  results  presented  above  suggest  that  a  continuation  of  the  perturbation 
solution  can  yield  a  more  accurate  representation  of  the  flow  field.  Some  ini¬ 
tial  attempts  in  this  direction  are  presented  in  this  section.  It  should  be 
emphasized  that  the  impact  of  the  initial  assumptions  is  likely  to  be  of  crucial 
importance  in  determining  the  validity  of  the  higher-order  correction  terms. 

In  particular,  neglect  of  viscous  terms  and  the  influence  of  turbulence  and  the 
turbulent  transition  process  in  the  formulation  may  make  such  efforts  of  aca¬ 
demic  interest  only.  It  is  the  opinion  of  the  writer  that  attempts  to  p.tduce 
a  theory  for  motors  with  large  L/D  should  incorporate  a  turbulent-  flow 

model  in  a  rational  way.  Methods  analogous  to  development  of  turbulent  pipe 
flow  might  be  appropriate.  Considerable  numerical  work  has  been  accomplished 
in  this  area,  notably  the  work  of  Beddini  (Refs.  13  and  14),  which  indicates 
the  importance  of  viscous  elements  of  the  flow  field.  Other  experimental  studies 
(Ref.  15)  also  suggest  a  similar  conclusion,  although  some  features  of  the  test 
procedure  and  their  influence  on  turbulent  properties  of  the  flow  are  in  ques¬ 
tion.  It  is  quite  likely  that  the  most  desirable  feature  of  the  present  analy¬ 
sis  -  its  simplicity  -  will  most  certainly  be  lost  if  viscous  effects,  especially 
turbulence,  are  explicitly  included.  In  spite  of  these  considerations,  the 

attempt  is  made  in  what  follows  to  extend  the  inviscid  model  to  0(M.  ).  It  will 

0 

be  shown  that  the  corrections  are  apparently  only  important  for  large  length  to 
diameter  ratio  combustion  chambers. 
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The  stream  function  to  O(M^)  Is  governed  by 
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and  the  momentum  balance  to 
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where  P  ^  -  ir2rV^  and  and  y^  are  given  by  Equations  A9  and 

A14  as  before.  The  boundary  conditions  are 
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A  solution  which  satisfies  these  requirements  is 
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The  corresponding  velocity  components  are 
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Thus  the  velocity  field  is  determined  to  0(M^  ).  Notice  that  only  the  axial 
correction  is  likely  to  be  of  any  consequence;  in  fact,  an  appropriate  approxi¬ 
mation  is 
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since  the  other  terms  are  of  order  zM^  or  smaller.  Thus  the  corrected  center- 
line  pressure  distribution  is  to  good  approximation: 
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in  physical  variables.  For  typical  chamber  parameters,  the  correction  is  of 
the  order  of  23  percent  for  (-^-)  =  kO. 


